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Critical hysteresis in ferromagnets is investigated through a iV-component spin model with ran- 
dom anisotropies, more prevalent experimentally than the random fields used in most theoretical 
studies. Metastability, and the tensorial nature of anisotropy, dictate its physics. Generically, ran- 
dom field Ising criticality occurs, but other universality classes exist. In particular, proximity to 
O(N) criticality may explain the discrepancy between experiment and earlier theories. The uniaxial 
anisotropy constant, which can be controlled in magnetostrictive materials by an applied stress, 
emerges as a natural tuning parameter. 



Hysteretic properties of ferromagnetic materials have 
long fueled applied research and, more recently, much 
theoretical interest Q- As a manifestation of the non- 
equilibrium dynamics of a disordered system with many 
degrees of freedom, hysteresis is described naturally in 
the language of statistical mechanics. A central aim of 
theoretical studies of hysteresis is to elucidate the ways 
in which microstructural details, such as domain config- 
urations, lattice structure, impurities or defects, affect 
macroscopic properties such as the shape of the hystere- 
sis loop and the Barkhausen noise statistics. 

A non-equilibrium version of the zero temperature ran- 
dom field Ising model (RFIM) has served to illustrate 
the competing effects of disorder and (ferromagnetic) ex- 
change interaction involved in hysteresis |3( • In three and 
higher dimensions, at weak disorder the model exhibits 
a discontinuous hysteresis loop, which becomes continu- 
ous at strong disorder. These two phases are separated 
by a critical loop for a given value of the disorder; as 
the latter is approached from the weak disorder side, the 
macroscopic discontinuity vanishes continuously, result- 
ing in a critical point characterized by universal scaling 
laws 0010 The corresponding critical exponents were 
obtained within a mean field approximation pertur- 
batively in a renormalizationgroup treatment and 
exactly on the Bethe lattice (While the model was 
originally suggested in part to relate this disorder induced 
critical scaling to Barkhausen noise measurements 0, Q , 
it seems that in most experiments the statistics of the 
noise is controlled instead by the depinning transition of 
domain walls 0, , which do not emerge simply from an 
analysis in terms of a RFIM.) 

Direct experimental evidence of disorder induced tran- 
sitions in ferromagnets was obtained only recently. A 
temperature controlled transition was reported for Co- 
CoO bilayers [j| and a similar transition was observed 
in Gd/W films subjected to different annealing proce- 
dures, which induce variations in the disorder through 
variations of the grain size [lflj. A study of hysteresis 
loops of Cu-Al-Mn alloys for different Mn concentrations 
and temperatures also identified a transition and 
the measured scaling exponents are consistent with those 
observed for Co-CoO bilayers. They do not agree, how- 



ever, with predictions of the RFIM. A natural explana- 
tion proposal for this discrepancy focuses on the nature of 
disorder; indeed, while random fields are convenient for 
theoretical exploration, they are seldom present in real 
ferromagnets, which display more complicated forms of 
disorder. Prominently, random anisotropies are present 
in most ferromagnets, including soft materials, and are 
believed to be particularly relevant in amorphous rare 
earth alloys (l2l |. 

A disorder induced phase transition was observed nu- 
merically in a random infinite anisotropy model, with ex- 
ponents close to the RFIM ones ^3|, supporting general 
symmetry arguments that were put forth in favor of uni- 
versality Q. However, infinite anisotropies pin the spins 
to given (random) directions, making each spin Ising-like 
on its own; as a result, the model is equivalent to a ran- 
dom field, random bond model. Furthermore, within a 
non-equilibrium context symmetry arguments ought to 
be taken with a grain of salt; it is known, for example, 
that the magnetization may point away from the applied 
field out of equilibrium, while in equilibrium minimiza- 
tion of the free energy requires alignment of the two. 
Such phenomena are a consequence of the presence of 
many metastable states (involved in the dynamics), and 
more systematic analyses that clarify their role and sub- 
stantiate the symmetry arguments are worthwhile. Along 
these lines, a renormalization group study of a random 
field vectorial (O(N)) model (RFVM), taking metasta- 
bility into account, showed that while one is justified 
in expecting a critical behavior identical to that of the 
RFIM generically, by tuning additional parameters dif- 
ferent universality classes may be visited |14| . 

Here, we analyze a non-equilibrium random anisotropy 
vectorial model (RAVM), in which iV-component spins 
are subjected to ferromagnetic interactions and random 
(finite) anisotropies. The N = 2 case was proposed in the 
past as a model of rare earth alloys, and its hysteretic 
behavior was studied numerically [l5|]- While in these 
studies the anisotropy averages to zero, here we allow 
for a fixed uniaxial component in addition to a random 
background. From the zero temperature spin dynamics, 
we construct the appropriate non-equilibrium effective 
action 0] which describes the evolution of the magnet 
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along the hysteresis loop. The most notable consequence 
of metastability is the generation of a 'random field term' 
in the action; in addition, the latter breaks the rotational 
symmetry verified by the equilibrium action. As a result, 
random anisotropy magnets indeed generically display 
usual, RFIM exponents (at least within the domain of 
validity of the perturbative analysis) for given values of 
the applied field and disorder strength. However, there 
exists a number of additional universality classes, and in 
particular a critical point with O(N) exponents which is 
reachable upon tuning of an additional parameter. The 
higher (tensorial) nature of the disorder in the RAVM 
provides such an additional parameter, namely the uni- 
axial anisotropy constant, in a natural fashion. As ex- 
plained below, proximity to a vectorial (O(N)) critical 
point may help explain the discrepancy between experi- 
mentally measured exponents and Ising ones. 

In the RAVM, A-component spins s on a d dimen- 
sional lattice interact via ferromagnetic nearest neighbor 
couplings Jo and a spin at site i couples to an anisotropy 
tensor K"^ '. In addition, the spins are subjected uni- 
formly to an applied magnetic field H, which varies (adi- 
abatically) in time and hence forces the system out of 
equilibrium, through a Hamiltonian 



H = - Jo ^2 ^ 
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where S{ • Ki • Si stands as a shorthand for 
Yla 8=1 s ?-^T^ s f • (Latin indices label lattice sites while 
Greek indices label spin components.) The anisotropy 
tensor may be decomposed into non-random and random 
components, and in the simplest (uniaxial) case 



Kf = K Q n c 
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where Kq is the uniaxial anisotropy constant and h a unit 
vector lying along the easy magnetization axis. The ran- 
dom components 8K a " are uncorrelated Gaussian ran- 
dom numbers with vanishing mean and standard devi- 
ation i?, so that the anisotropy tensors are distributed 
according to the density 
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The parameter R plays an analogous role here to that 
of the width of the disorder distribution (also called R) 
in random field models 0, Q|. Finally, for the sake 
of calculational simplicity, instead of fixed length spins 
(with, e.g., \si\ 2 = 1 for each site i) we consider 'soft 
spins' whose lengths can take any values. Following Refs. 
PL 13 ■ f° r stability we add to the Hamiltonian a sum 
of single site terms ^2iV(si), so that, at each site, a 
Mexican hat potential V(si) — —a \si\ 2 /2 + b |s*j| 4 /4 pre- 
vents the spin from diverging. In the appropriate limit 
(a = b — > oo), soft spins reduce back to unit spins, but 
as the length of spins is modified under renormalization, 
the specific (bare) values of a and b are irrelevant. 



As mentioned, the applied field H varies (adiabat- 
ically) in time and forces the spins through a non- 
equilibrium trajectory. In order to study the critical be- 
havior of the system, we may confine ourselves to the 
simple zero temperature relaxational dynamics 
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where V is an effective damping coefficient. We point out, 
though, that Eq. @ is certainly not the most realistic 
choice of dynamics, which in general includes precession 
of spins and a more complicated damping factor, better 
described by a Landau-Lifshitz-Gilbert-like equation @|. 
Nevertheless, we expect the critical behavior not to bear 
crucially on the specifics of the microscopic dynamics, 
and Eq. IgJ appears as the simplest candidate for an an- 
alytic treatment. In the same vein, in what follows we 
make a final, customary simplification in replacing the 
lattice spins Si by a continuum vector field s(x). 
In the continuum approximation, Eq. Q becomes 



r f?ijM = JV 2 s + H(t)+K( X ) 



s + as — b\s\ s, (5) 



where the constant J results from the continuum ex- 
pansion of the exchange interaction 0] , K (x) • s 
is a shorthand for the vector field with components 
Y^3—i ( x ) s<3 j and higher orders (in derivatives and 
possibly fields) have been neglected. Aspects of the 
hysteretic critical behavior are more transparent in the 
language of a generating functional 16] than directly 
through the equation of motion. In the usual fashion 
0, Il4j , we introduce an auxiliary field <fi to exponenti- 
ate a delta function that forbids any trajectory that does 
not obey the equation of motion, resulting, up to some 
constant prefactors, in a functional 



X>sX>0exp ( / dtd d x 



,9s 
di 



+JV 2 s + H + K-s + as-b\sfsjj (6) 

which captures the possible histories of the system. The 
advantage of this procedure is that now one can easily 
average the generating functional over the distribution 
p(K a P) for anisotropy tensors K a @ at all positions, as 



VK a Pp(K al3 )Z 

VsVcj) exp (S'cff [s(x,i),0(x,i)]), (7) 



with an effective action 



JV s + H + K-s 



S'cff = fdtd d x{c, 
+as-b\sf s + K n(ii-s) + ^ J dt'(<p '■ $')(s- s*)} 
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where s, cf> and s",</>' are evaluated at times t and t', 
respectively. The effective action encodes the averaged 
solutions of Eq. © and avoids one the complication 
of solving a stochastic equation first and then averag- 
ing. In carrying out the average, one trades the stochas- 
tic (anisotropy) term with new terms coupling s and <fi, 
and the quenched nature of the disorder is reflected in 
the presence of a double integral over time (over t and 
t'). In the RAVM, the two terms with (<fi ■ n) (n ■ s) and 
(cj)-(f>')(s-s') replace random field terms of the form (</>•</>') 
in the RFIM || and the RFVM 

As mentioned, the effective action in Eq. (JHJ encom- 
passes all the solutions and as a result is invariant under 
the transformation (H,s,(f)) — ► (— if, — s, — <^>). The hys- 
teresis curve, however, is not symmetric in general upon 
inversion of the magnetic field and magnetization; in par- 
ticular, the value of the magnetization at zero field (re- 
manent magnetization), and vice versa that of the field 
when the magnetization is zero (coercive field), do not 
vanish. This is because the system follows in reality 
a given metastable state which evolves along with H(t) 
(more precisely, each branch of the hysteresis loop corre- 
sponds to a metastable trajectory). Following a trick of 
Ref. Q , we get rid of the unwanted solutions by shifting 
the field s by a quantity <?(i), which represents the av- 
eraged sum of all the 'unphysical minima'; the resulting 
effective action, S s, metastable, encapsulates the magneti- 
zation and response function along a branch of the hys- 
teresis loop. For the sake of simplicity, we consider first 
the case in which the field H is applied along the easy 
axis of magnetization, i.e., H — Hh; we comment below 
on the general case. We then expect the magnetization, 
and hence the vector <r, to lie along h too. Shifting the 
spin field according to s — *■ s + a(t)h, we obtain 



-r 



3s|| 

m 



ScS, metastable = / dtd d X j(/>|| 

+ JV 2 S|| + H + (s|| + a) (a-b \s 

■ [-rfy + JV 2 

■Jdtdt'd d x U\\</>\\ 
\s {] +a(t)) (*f | +<7(f)) + SW_ L 



r 



01 



III 

2 



s± + s±ia 



)]} 



i 9 ) 

where we have decomposed s = (sii , s±) and <f> = (^>ii , 
into longitudinal and transverse components with respect 
to the direction given by ii. 

A number of results may be deduced from the form of 
the corrected action S e ff, metastable- In Eq. ©, the bare 
longitudinal 'mass' a (the coefficient of the </>||S|| term) is 
dressed into aii = a — 3ba 2 + Kq and the bare field H 
into H = H — Tdt<J + a(a — ba 2 + Kq). As, in general, 
ay and H do not become small (or vanish) simultane- 
ously at H = 0, criticality does not occur at vanishing 
field. This reflects the non-equilibrium nature of the tra- 
jectory, chosen among many metastable states generated 
by the disorder. The more remarkable manifestation of 



metastability is, however, the generation of an effective 
'random field' a(t)a(t')cj) • 4>' term. As a result, the ef- 
fective action in Eq. JHJ differs from its random field 
counterpart by the presence of additional terms of the 
form (j>cf)'s^ and <p(j)' ss' . Without these terms the ac- 
tion exhibits a non-trivial (non-Gaussian) critical point 
below the upper critical dimension d c = 6, which can 
be characterized by a perturbative renormalization group 
treatment in d = 6 — e dimensions Now, power 

counting predicts that this critical point is stable with 
respect to the extra random anisotropy terms: from the 
natural rescalings x — > bx, t — > b 2 t, s — > b 2 ~ d / 2 s, and 
4> — > 6~ 2 ~ d / 2 0, we find a scaling dimension of (4— d) /2 for 
the (jxfi's^ terms and of 4 — d for the cjxjj'ss' terms, which 
are thus irrelevant close to d c = 6 dimensions. Conse- 
quently, at least within the perturbative domain, critical- 
ity in the RAVM is identical to that in the RFVM 
with, generically, RFIM exponents |4| reflecting 'mass- 
less' fluctuations of the longitudinal component sii. 

By symmetry, as in the RFVM there must exist 
here aO(JV-l) critical point representing 'massless' fluc- 
tuations of the transverse components s± , corresponding 
to spontaneous magnetization in the transverse direction. 
Then, upon appropriate tuning of an additional param- 
eter, these two critical points may merge, resulting in 
a rotationally invariant vectorial (O(N)) critical point. 
The latter occurs for a symmetric action, in particular 
the effective longitudinal mass a|| and the effective trans- 
verse mass a± = a — ba 2 must become small simulta- 
neously. In the RFVM, simultaneous vanishing of the 
effective masses and applied field is possible only if the 
magnetization vanishes at H = 0, i.e., for very 'thin' hys- 
teresis loops (with small area) Here, crucially, the 
higher (tensorial) nature of anisotropy alters this picture: 
vector criticality occurs generically at non-vanishing val- 
ues of the applied field and magnetization. Indeed, since 
Kq modifies both the field and the longitudinal mass, the 
values of H, a\\, and a± may become critical simultane- 
ously at non- vanishing values of H and a. As a result, 
in the RAVM the hysteresis loop need not be 'thin' to 
display vectorial criticality. The uniaxial anisotropy con- 
stant Kq may be tuned instead of the disorder width to 
reach Ising criticality, or along with the disorder width 
to reach vectorial criticality. 

We emphasize that this picture is a direct consequence 
of metastability. By contrast, the equilibrium random 
anisotropy model |17| displays a lower critical dimension 
of d c = 4, which is also, if naive power counting is to be 
believed, the upper critical dimension. Systematic stud- 
ies of the behavior about d c = 4 are plagued with a num- 
ber of technical difficulties and, although the weak 
disorder phase below d c was captured in a recent analyt- 
ical treatment |19| , agreement with experiments |2(| and 
numerics [2l| is still controversial. While it is certainly 
legitimate to ask whether similar difficulties arise out of 
equilibrium away from d c = 6, we note that, curiously, at 
least in the perturbative domain, metastability simplifies 
the problem. 
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So far, we have considered the particular case in which 
the system is magnetized along its easy axis. A similar 
analysis may be applied to the general case, in which the 
magnetization lies along a direction \x intermediate be- 
tween those of n and H . As H increases (or is varied) in 
time, jl rotates in space. This, however, does not affect 
the analysis significantly (calculationally, because differ- 
ent times effectively decouple in the action). The fact 
that jl does not lie along n or H changes the symme- 
tries, in particular because additional transverse terms 
are generated by the shift in a. These seem to allow 
for the possibility of 0(2) and 0(N — 2) critical points, 
provided additional parameters may be tuned. 

In sum, we have shown that critical hysteresis in the 
RAVM is described, generically, by RFIM exponents. 
Thus, we expect the conclusions of simulation studies of 
the random infinite anisotropy model to extend to 
the finite anisotropy case in general. However, we also 
expect it to be easier to identify vectorial critical points 
in the presence of anisotropy than in the presence of ran- 
dom fields, and a potential explanation for the discrep- 
ancy between experiments and theory lies in a putative 
proximity of the regime in which experiments are carried 
out to such a vectorial critical point. An experimental 
study in which various parameters are scanned systemat- 
ically should reveal whether proximity to a vectorial crit- 
ical point is verified. If it is the case, one expects O(N) 
exponents, that cross over to RFIM ones only above some 
scale which may be rather large 0, ^| . 

For practical reasons, the standard experimental tun- 
ing parameter is temperature, either annealing or mea- 
sured temperature, but a number of interpretation prob- 
lems are associated with the corresponding techniques. 
In the first case, one loses control over changes in mi- 



crostructure and one is bound to repeat the experi- 
ment on different samples, with possibly large sample-to- 
sample variations. In the second case, varying the tem- 
perature modifies simultaneously several physical quanti- 
ties. In the RAVM, the natural tuning parameter is the 
uniaxial anisotropy constant Kq, rather than the tem- 
perature; it may be used both to reach the RFIM and, 
along with a second tuning parameter, to look for vec- 
torial criticality. Tuning Kq seems a good experimental 
possibility for magnetostrictive materials in which, in the 
simplest description, an applied stress r along the easy 
axis shifts the value of the uniaxial anisotropy from Kq 
to Kq +3At/2, where A is the magnetostriction constant 


A more drastic reason for the discrepancy between 
experiment and theory may be, of course, that the ex- 
perimentally predominant form of disorder is neither of 
the random field nor of the random anisotropy type. In 
particular, in materials used at present in experiments 
H> El E3 , the presence of random bonds and demagne- 
tizing fields might alter the theoreticalpicture Q, and 
so would putative strong dipolar forces [Tj. Nevertheless, 
we expect our results to be relevant for a wide class of 
amorphous ferromagnets, and in particular it would be 
interesting to check them against experiments on poly- 
crystals, in which dipolar forces are weak and anisotropy 
is the dominant form of disorder. 
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